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Th 1 ' ABSTRACT 



o 



This paper explores the quantitative connection between globular clusters and the "dif- 
fuse" stellar population of the galaxies they are associated with. Both NGC 1399 and NGC 
4486 (M87) are well suited for this kind of analysis due to their large globular cluster popula- 
tions. 

The main assumption of our Monte Carlo based models is that each globular cluster 
is formed along with a given diffuse stellar mass that shares the same spatial distribution, 
chemical composition and age. The main globular clusters subpopulations, that determine 
the observed bimodal colour distribution, are decomposed avoiding a priori parametric (e.g. 
Gaussian) fits and using a new colour (C-Ti)-metallicity relation. The eventual detectability 
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\ of a "blue" tilt in the colour magnitude diagrams of the blue globulars subpopulation is also 

• addressed. 

A successful link between globular clusters and the stellar galaxy halo is established by 

assuming that the number of globular clusters per associated diffuse stellar mass t is a function 
of total abundance [Z/H] and behaves as t — 7 exp(— S[Z/H] ) (i.e. increases when abundance 
decreases). 

The simulations allow the prediction of a surface brightness profile for each galaxy 
through this two free parameters approximation. The 7, 6 parameters that provide the best 
k>( \ fit to the observed profiles in the B band, in turn, determine several features, namely, large 

• scale halo colour gradients, globular clusters-halo colour offset, clusters cumulative specific 

frequencies, and stellar metallicity distributions, that compare well with observations. 

The results suggest the coexistence of two distinct stellar populations characterised by 
widely different metallicities and spatial distributions. One of these populations (connected 
with the blue globulars) is metal poor, highly homogeneous, exhibits an extended spatial dis- 
tribution and becomes more evident at large galactocentric radius contributing with some 20% 
of the total stellar mass. In turn, the stellar population associated with the red globulars is ex- 
tremely heterogeneous and dominates the inner region of both galaxies. 

Remarkably, and although the cluster populations of these galaxies exhibit detectable 
differences in colour distribution, the 6 parameter that determines the shape of the brightness 
profiles of both galaxies has the same value, 6 « 1.1 to 1.2 ±0.1. 

Key words: galaxies: star clusters: general - galaxies: globular clusters: individual: NGC 
1399, NGC 4486 - galaxies:haloes 
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1 INTRODUCTION 

The idea that globular clusters (GCs) harbour important clues in re- 
lation with the early stages of galaxy formation is a widely accepted 
concept. One of the most compelling arguments in favour of the 
existence of a connection between GCs and major star formation 
episodes in the life of a galaxy is the constant cl uster formation effi - 
ciency, defined in terms of total baryonic mass dMcLaughlitJ 19991) . 
in different galaxies. 

However, breaking the code that leads to a detailed quan- 
titative link between GCs and the underlying "diffuse" stel- 
lar population is still an open questio n. Such a connect ion 
has been disc ussed on theoretic al (e.g. iBeaslev et al.l |2002| or, 
more recently, Pipino et"al] |2007) and observational grounds (e.g. 



Forbes & Forte 



200lh. In t he particular case of Milky Way GCs, 



Pritzl. Venn & Irwirj J2005I) ~ound that the chemical similarities be- 
tween clusters and field stars with [Fe/H] sC — 1 suggests a shared 
chemical history in a well mixed early Galaxy. 

Clarifying this issue may certainly yield some arguments in 
favour (or against) some predom inant ideas that have been widely 
refere nced in the literature (e.g. Egg en et "alii 19621 : ISearle & Zinnl 
Il97£ ) and later explored w it hin the frame of d ifferent scenarios 
(e.g jAshman & Zepj|l992l or lForbes etaljl997h . 

A good perspective of the com plex situation in t his c ontext 
is given in th e thorough review bv lBrodie&Stradeij lEofji) and 
iLarsonl d2006l) . 

An initial confrontation between GCs and halo stellar popula- 
tions shows more differences than similarities: a) In general, GCs 
exhibit more shal low spatial distribu t ions than those cha racterising 
galaxy light (e.g. lRacine et al.ll 19781: iDirsch et alj|2003h : b) There 
is a colour offset in the sense that mean integrated globular colours 
appear bluer tha n those of the gal a xy halos at the same galac- 
tocentric radius 1 Str om et al. 1 1 19811 : iForte. Strom & Stroml fl98ll : 
Ijordan et alj|2004 ): c) GCs show frequently bimodal colour (and 
hence , metallicity) distributions (see, for example, IPeng et all 
l2005t) . This feature does not seem exactly shared by the resolved 
stella r populations in nearby resolved galaxies (se e iDurrell et al 



2001 



2006) 



lHarris & Harris! 120021 : iReikuba et all 120051 or 



Mouhcinc 



As discussed later in this work, those differences arise, 
mainly, from the fact that GCs analysis usually provide number 
weighted statistics while galaxy halos observations yield luminos- 
ity weighted measurements. 

A preliminary quantitative ap proach to the globulars - stellar 
halo connection was presented in IForte, Faifer & Geislen (2005) 
(hereafter FFG05). This last paper shows that, given the areal den- 
sity distribution of the "blue" and "red" globular cluster subpopula- 
tions in NGC 1339, the galaxy surface brightness profile, galacto- 
centric colour gradient and cumulative GCs specific frequency, can 
be matched by linearly weighting the areal density profiles. The 
"weight" corresponding to each component of the brightness pro- 
file is the inverse of the intrinsic GCs frequency characteristic of 
each cluster population. 

The main argument behind that approach is that the shape of 
the colour (and metallicity) distribution of each globular cluster 
subpopulation does n ot change with galactocentric radius. La rge 
angular scale studies jDirsch et alj 120031 : iBassino et al.l 120061) in 
fact show that the colour peaks in the GCs colour statistics of NGC 
1399 keep the same position (or show very little variation) over 
large ga lactocentric ranges. A similar result is obtained for NGC 
4486 by Kundu & Zepf ( 2007) who find those peaks do not show a 
detectable variation in colour over 75 kpc in galactocentric radius. 



It must be stressed that those subpopulations are "phenomeno- 
logically" defined in terms of their integrated colours but each 
might eventually have a given spread in age and/or metallicity. 

The presence of a "valley" in the globular colour statistics is 
usually adopted as a discriminating boundary between both sub- 
populations. The need to revise such a procedure was already sug- 
gested by figure 5 in FFG05. This last diagram showed that the 
NGC 1399 GCs bluer than the blue peak have a distinct behaviour 
of the areal density profile, exhibiting a flat core that disappears 
when all "blue" clusters (i.e. all GCs bluer than the colour valley) 
are included in the sample. That result prompted for a further anal- 
ysis as discussed below. 

This paper generalises the FFG05 approach trough Monte- 
Carlo based models. In this frame, "seed" globulars are generated 
following a given abundance Z distribution and then associated 
with a "diffuse" stellar mass that shares its age, chemical compo- 
sition and statistical spatial distribution. The luminosity associated 
with this mass is derived from a mass to luminosity ratio adequate 
for a given age and metallicity. These models aim at reproducing 
the features mentioned above and seek for a function that could link 
GCs and diffuse stellar populations keeping a minimum number of 
free parameters. 

Both NGC 1399 and NGC 4486 appear as adequate targets in 
order to perform such a modelling due to their prominent globular 
cluster systems (GCS). Although these systems show some struc- 
tural similarities in terms of their spatial distribution, they also dif- 
fer markedly both i n the shape of the ir GCs colour statistics and 
specific frequencies dForte et alj2002b . 

This work also presents new Washington photometry, ob- 
tained and handled in a homogeneous way, that allows for a re- 
discussion of the GCS properties in the inner region of both galax- 
ies and, in particular, of the behaviour of the areal density of GCs 
with colour. In turn, recent wide field photometric stu dies of the 



GCS associated with NGC 1399 dBassino et al.Jl2006l) and NGC 
4486 dTamura et al]|2006d : iTamura et alj|2006bh , are well suited 



for extending the analysis to larger galactocentric radii. 



2 OBSERVATIONS AND DATA HANDLING 

Photometric observations were carried out with the Mayall and 
Blanco 4-m telescopes at KPNO and CTIO respectively and 2048 
pixels on a side CCDs, with a pixel scale of 0.43 arcsecs. T he C 
filter from the Washington System dHarris & Canternal ll979l) and 
the Rkc filter of the Kron-Cousins system were used at both tele- 
scopes. As noted by iGeislerl d 19961) this last filter is comparable 
to the Ti filter in the Washington system although much more ef- 
ficient in terms of transmission. In what follows we keep the Ti 
denomination for our red magnitudes. 

Two Rkc images (exposure: 600 sees each) and three C im- 
ages (exposure: 1500 sees each) were secured for both galaxy 
fields. 

Seeing quality for these frames varied between 1.0 and 1.6 
arcsec. These images were processed with the CCDRED routines 
within IRAF0 including bias and flat fielding. 

The galaxy background removal was performed with routins 
included in the VISTA image prossesing systems. 



1 IRAF is distributed by the National Optical Astronomical Observato- 
ries, which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation 
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Figure 1. Distribution GCs candidates brighter than Ti = 23.2. Left panel: NGC 1399 field. North is to the right and East is up. Right panel: NGC 4486 field. 
North is up and East to the left. In both panels, circles have radii of 120 and 420 arcsec respectively. The GC sample completeness in these areas is close to 
95%. 



Table 1. Photometric data NGC 1399. 



ID X (arcsec) Y (arcsec) Ti (C-Ti) roundness 



537. 196.8 -427.6 24.15 0.63 0.77 

553. 136.5 -427.5 24.33 0.65 0.85 



Table 2. Photometric data NGC 4486. 



ID X (arcsec) Y (arcsec) Ti (C-Ti) roundness 



23314. -13.1 -126.5 21.67 1.43 0.93 

23406. -32.1 -124.9 23.39 1.26 0.94 



PSF photometry on all the frames was then car ried out using 
the ALLFRA ME version of the DAOPHOT package dStetsonl 1987k 
IStetsodll99lh . 

The instrumental photometry was trans formed to the s tandard 
system using calibration standard stars from lGeislej {l996). 

Image classification, in terms of resolved and non-resolved 
objects, was performed as in iForte et alj d200ll) . Briefly, that pro- 
cedure combined the use of the round and sharpness parameters 
defined in DAOPHOT and also the mirrored envelope of the Ti 
PSF ALLFRAME magnitude vs. the difference between this mag- 
nitude and that obtained using aperture photometry for every object 
detected on the images. 

Non-resolved objects brighter than Ti = 23.2 and with (C- 
Ti) colours between 0.9 and 2.3 were considered as cluster candi- 
dates and their distribution on the sky is depicted in Figure[T] 

Circles with r=120 and 420 arcsec delineate the area used for 
the analysis of the surface density distributions in the inner regions 
of the galaxies. 

Figure [2] shows the errors on the (C-Ti) colours as a function 
of Ti magnitude as delivered by DAOPHOT. A median error for 
the (C-Ti) colours of ± 0.07 mags, is reached at Ti = 23.2, which 
is adopted in what follows as the limiting magnitude of the analysis 
in order to assure good quality colours. 

The photometric data for both galaxies, Table Q] and [2] are 
available in the electronic journal version. Coordinates are referred 
to the galaxy centers and defined as in Figure Q] 

ADDSTARS experiments were carried out to estimate the 
completeness of the non-resolved objects (which is expected to be 
the case for GCs at the distances of NGC 1399 and NGC 4486 from 
the sun). Ten images, adding 1000 artificial stars each, on both C 
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Figure 2. (C-Ti) colour errors as a function of Tj magnitude for both 
galaxy fields. The vertical line at Ti = 23.2 is the limiting magnitude 
adopted in the analysis. The median colour error for the sample is 0.04 
mags. 



and Ti master images, yielded a completeness of 94 and 96% at 
Ti = 23.2, for NGC 1399 and NGC 4486 respectively. 

A comparis on field of 77.7 arcmin was taken from 
Forte et al I diooih who performed C and Ti photometry following 
the same procedure. This field has 146 non-resolved objetcs within 
the colour-magnitude boundaries adopted for the globular cluster 
candidates. 



3 COLOUR-MAGNITUDE DIAGRAMS AND COLOUR 
DISTRIBUTIONS 

The Ti vs. (C-Ti) colour diagrams for non-resolved objects are 
displayed in Figure|3] The limiting magnitude Ti = 23.2 is indi- 
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Figure 3. Ti vs. (C-Ti) colour magnitude diagram for unresolved objects 
in the field of NGC 1399 (upper panel) and NGC 4486 (lower panel). The 
domain of the GCs candidates discussed in the text is indicated by the rect- 
angular area. The tilted line, for NGC 4486, is defined by the modal values 
of the colour statistic inside 0.5 mags, intervals in Ti for the blue GCs. 



cated as a horizontal line while vertical lines at (C-Ti)=0.90 and 
(C-Ti)=2.30 define the domain of the globular cluster candidates. 

A distinctive feature in the lower panel of Figure[3] in contrast 
with the upper panel of the same figure, is a noticeable tilt of the 
colours of the blue clusters associated with NGC 4486, which is not 
detectable in the case of NGC 1399. The tilted line Figure [3] (lower 
panel) was obtained by fitting the modal values on a smoothed im- 
age of the colour magnitude diagram (adopting a round Gaussian 
kernel of 0.05 mags.) yielding: 



Ti = 41.50 - 16.67(C - Ti) 



(1) 



This relation implies a 0.06 mags. (C-Ti) colour increase per 
magnitude that is compar able to that detectab le in the (g-z) colours 
of 0.045 per magnitude dStrader et al. 2006). The possible reason 
for the presence of a "blue" tilt in the colour magnitude diagram is 
discussed in Section[8] 

The (C-Ti) colour histograms for GCs within a circular galac- 
tocentric region defined between 120 and 360 arcsec are depicted 
in Figure [4] These histograms have been corrected for back- 
ground contamination by subtracting the comparison field his- 
togram (scaled by area and also shown in these figures) and contain 
about 1000 and 1800 candidate GCs for both galaxies, respectively. 

In order to minimise the presence of bright objects tha t 
might be identified as compact dwarfs (see IPhillips et alj 1200 ll) 
we adopted an upper cut off at Ti=21.0. We stress that 
lOstrov. Forte & Geislen 1 d 19981) noted that NGC 1399 cluster can- 
didates brighter than this magnitude exhibit a unimodal colour dis- 
tribution, a feature later confirmed in lDirsch et al.1 ( 120031) . As a ref- 
erence we point out that Omega Cen-like objects would appear at 
Ti « 21.4 and 21.0 for NGC 1399 and NGC 4486, respectively. 

These last figures also indicate the position of the so called 
"colour valleys" at (C-Ti)=1.55 and 1.52 for NGC 1399 and NGC 
4486 respectively. These values were determined using Gaussian 
smoothed colour histograms with a colour kernel of 0.05 mags. The 
same procedure leads to (C-Ti)=1.26 and 1.21 for the blue peaks 
and (C-Ti)= 1.75 and 1.72 for the red peaks in both galaxies. 



4 DESCRIPTION OF THE MODEL 

In this section we describe each of the steps involved in the model 
an the main hypothesis behind it, namely: 

a) The decomposition of the colour histograms in terms of the 
cluster subpopulations leading to their [Z/H], [Fe/H] and (C-Ti) 
colour distributions. 

b) The determination of the projected areal density distribution 
for each of the cluster subpopulations. 

c) Establishing the link between each cluster and its associated 
diffuse stellar population. 

d) Deriving the parameters that determine the shape of the pre- 
dicted galaxy surface brightness profile. 

a) The decomposition of the colour histograms. 

The first step is the decomposition of the the observed colour 
histograms shown in Figure [4] avoiding an a priori functional de- 
pendence (e.g., the usual Gaussian assumption). It must be empha- 
sised that, matching the two-peaked colour histograms observed 
both in NGC 1399 and NGC 4486 through the adopted colour- 
metallicity relation (see below), necessarily requires two distinct 
globular cluster populations. "Seed" clusters were then randomly 
generated in the abundance Z domain according to a given statisti- 
cal dependence. Trial and error shows that exponential behaviours 
f(Z) w exp[— (Z — Zi)/Zs] where Zs is the abundance scale 
length and Zi is the minimum abundance, provides acceptable fits 
to the observed histograms (i.e., within the Poissonian uncertainties 
associated with each statistical bin). Some more complex functions 
cannot be rejected but would imply a larger number of free param- 
eters not justified in terms of those uncertainties. 

As for the minimum abundance, we adopted Zi = 0.003Zq , 
that corresponds to [Fe/H] = —2.65 in the metallicity scale pre- 
sented below. As discussed later, a dependence of Zi with Ti leads 
to a blue tilt that reproduces the colour-magnitude diagram of the 
NGC 4486 GCS. 
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Figure 4. (C-Ti ) background corrected colour histogram for NGC 1 399 (left panel) and NGC 4486 (right panel) globular candidates within a circular annulus 
defined between 120 and 360 arcsec in galactocentric radius. In both cases the areal scaled subtracted background is shown by the dotted line histograms. 
The combined counting statistical error bars are also shown. The histograms contain fa 1000 and fa 1800 cluster candidates for NGC 1399 and NGC 4486, 
respectively. Vertical lines idicate the position of the colour valleys. 



The decomposition procedure aims at matching the position of 
the colour peaks and colour valley while keeping a minimum valu e 
of the quality index of the fit, \ 2 , defined as in lCote et alj lfl998). 

The cluster abundance Z was linked to metallicity on the Zinn 
and West (1984) scale. The ado ption of the [Fe/H] z w scale imply 
some caveats (see for example, iThomas et all2003l or lStrader et al.l 

2007) about the nature of this index. In this work we use the rela- 

I 1 J — — ^J. 

tion found by Mendel, Proctor & Forbes (2007) for the stellar pop- 
ulation models given bv lThomas. Maraston & Ko rn (2004) : 

[Fe/H]zw = [Z/H] - 0.131 (2) 

An integrated colour was then obtained for each cluster 
through an empirical (C-Ti)-[Fe/H] colour-metallicity calibration. 

Several approach have been made in the past aiming at de- 
termining the colour met allicity relation. For exa mple, the original 
linear relation found by iGeisler & Forte] (T9 90) for MW clusters 
was later improved by [Harris & Harris! j2002). Being an important 
step in the modelling process, we attempted a new calibration, de- 
scribed in Section|5] which yields a quadratic relationship between 
metallicity and integrated cluster colours. 

Before comparing the model cluster colours with the observed 
histograms, we added interstellar reddenings, (E(B — V) = 0.015 
and 0.022 for NGC 1399 and NGC 4486, respectively) from the 
ISchlegel et al] dl998h maps, adopting E(C - Ti) = 1. 97 E(B - 
V), and Gaussian errors matching their behaviour as a function of 
cluster brightness displayed in Figure|2] 

Examples of the decomposition procedure are shown in fig- 
ures [5]and|6] The first dia gram displays the results obtained from 
fitting the kernel colour distribution of GCs belonging to the bright- 
est galaxies sample in the Virgo ACS jPeng et al]|2005l ; figure 5). 
These galaxies are comparable in brightness to both NGC 1399 
and NGC 4486. The parameters of the best fit are Zs(blue) = 
0.035Zq for the blue clusters (23% of the total population) and 
Zs(red) — 1.05Zq for the red clusters. In this case, we trans- 
formed the (C-Ti) colours to (g-z) by adopting: 

{g-z) = (C-Tr) -0.29 (3) 



600 
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Figure 5. Model fit (continuous line) to the kernel histogram representa- 
tive of globular clusters associated with the brightest ellipticals (black dots) 
in the Virgo ACS bv lPeng et al . (2005). The adopted Gaussian colour ker- 
nel is 0.05 mags. The two components correspond to Zs of 0.035Zq and 
1.05Zq for the blue (dotted line) and red (dashed line) clusters subpopula- 
tions respectively with peaks at (g-z)=0.92 and 1.43 (vertical lines). 



This colour transformation is consistent with the colour of the 
peaks in the ACS bright galaxies sample compared to our estimate 
of the (C-Ti) pe aks in Fig ure [4] and also with the colours relation 
derived from the Maraston (2004) models. 

Figure|6]shows th e best fit obtained for 91 M W globulars with 
(C-Ti) colours (from [Harris & Canternal [l977l) or (B-I) colours 
(transformed according to (C - Ti) = 1.03(S - I) Q - 0.43) 
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Figure 6. Upper panel: Colour histogram for 91 MW globular clusters. The 
bars indicate the statistical count uncertainties. The continuous line is a 
model fit described in text. Lower panel: Model components for MW glob- 
ulars. Blue clusters are fit with an abundance scale Zs(blue) = O.O35Z0 
and red globulars with Z sired) = 0.50Zq. Colour peaks at (C — Ti) = 
1.15 and 1.51 are indicated by vertical bars. 



in iReed et all dl988h . In this case, 70% are assigned to the blue 
clusters subpopulation, with Zs(blue) = 0.035 iO.OlZ©, and the 
remaining 30% to the red one, with Z sired) = 0.50 ± 0.05^0. 
These parameters imply [Fe/H] peaks at -1.7 and -0.5, respec- 
tively. The small sample of MW clusters has large statistical 
uncertainties but the fit is consistent with the observed shape of 
the [Fe/H] distribution (see lBica et alj20 06 and references therein). 



The model assumes that each of the GCs subpopulations 
has its own and distinctive spatial distribution. As noted before, 
however, the adoption of a given colour window to define each 
cluster subpopulation, leads to ambiguous results in the case of 
NGC 1399. This particular aspect is discussed with more detail in 
Section [6] on the basis of the photometry presented in this paper. 
In particular, we stress that the variation of the slope observed for 
the blue GCs (depending on the colour window adopted as their 
domain) may arise as a partial superposition of the two cluster 
subpopulations. 

c) T he GCs-diffuse stellar p opulation link. 

IZepf & Ashman! (fl99 3) introduced the T parameter, defined as 
the total number of GCs per galaxy stellar mass unit. In this work, 
we generalise that parameter by assuming that the number of globu- 
lar clusters per associated diffuse stellar mass t is a function of total 
abundance [Z/H]. After exploring different possible functions, we 
adopted: t — 7 exp(— 8[Z/H]) (i.e. t increases when abundance 
decreases), and then: 



dN/d[Z/H] = t([Z/H])M{[Z/H]) 



(4) 



where dN is the number of CGs associated with a stellar mass M 
and an abundance [Z/H] that belongs to a given subpopulation. This 
assumption leads to a "diffuse" stellar mass per cluster with a given 
[Z/H]: 



M* = 1/t 

and then to an integrated luminosity: 

L = M*/(M/L) 



(5) 



(6) 



Where (M/L) is the mass-luminosity ratio characteristic for stars 
with the same age and metallicity of the "seed" globular cluster. 

In this w ork we adopted the (M/L) ratio for the B (Johnson) 
band given bv lWorthevl < ll994l) and an age of 12 Gy. This ratio can 
be approximated as 



(M/L) B = 3.71 + ([Z/H] + 2.0) 2 



(7) 



This approximation differs from the Worthey's (M/L) ratios 
by, at most, « 7%. 

Note that we adopt [Z/H] instead of [Fe/H] as Worthey's mod- 
els assume solar scaled metallicities. A comparison with other 
models gives an idea about the u ncertainty in this ratio. For ex- 
ample, models in Maraston (2004), for the same age and a Salpeter 
luminosity function, show an overall agreement better than 10% 
with Worthey's except at the lowest abundance where they deliver 
a ratio ~ 24% larger. The effect of age variations on this ratio is 
discussed below. 

In particular, we choose the B and R bands since large scale 
surface photometry is available for both NGC 1399 and NGC 4486 
(see Sections |7] and [8) and no comparable data has been published 
in the C and Ti bands for both galaxies. 

Although th e (M/L)b ratio d epends on age, we stress that 
most works (e.g. Jordan et alj|2002h have not detected significant 
age diff erences for the clu ster subpopulati ons in NGC 4486 . 
In turn, Forbes et a D j200ll) (and also see IPierce et al] l2006al: 
iPierce et al] l2006bl or THempel et al] l2007h find arguments to 
support the presence of a fraction of "intermediate age" clusters in 
NGC 1399 and in other galaxies. However, age differences as large 
as ± 2 Gy will not have a strong impact on the integrated colours. 



b) Projected spatial distributions. 



d) The shape and colour of the galaxy surface brightness profile. 
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Each stellar mass element associated with a given "seed" GC 
(and determined by the adopted 7 and S parameters) was split into 
a number of "luminous" particles (i.e. 100 per cluster). These par- 
ticles were statistically located on the plane of the sky by adopting 
the same spatial distribution that characterises each of the cluster 
subpopulations in order to construct bi-dimensional blue image (2 
arcsec per pixel) of the galaxies. A red image was also obtained by 
transforming the (C-Ti) colour of each luminous particle to (B-R) 
by means of: 

(B-R)kc = 0-704(C-Ti)+ 0.269 (8) 

empir ically obtained f rom MW GCs with Joh nson dReed et all 
Il988h and Washington {h arris & Canterna|[l977h photometry. 

The synthetic B and Rkc images were then analysed with 
the task ELLIPSE within IRAF in order to derive surface bright- 
ness profiles and colour gradients along the semi-major axis of the 
galaxies and, in the case of NGC 4486, the variation of ellipticity 
along the same axis. 

This treatment generalises the profile expression given in 
FFG05: 

Ms = (V -Mv) + A(B)+2Mog\s B (red)\ 



',.5log 



+ 



N(blue) 
C B 



] 



(9) 



where N(blue), N(red) are the areal densities of the blue and 
red clusters at a given galactocentric distance and Cb = 
S (blue) /S (red). Introducing the definition of the t parameter 
given before, leads to: 



1.0 



[Z/HU 



dn 



1.0 



7 / m 4Z/H] (MjL) 



exp{S[Z/H]}d[Z/H] (10) 



where the integrals are performed on the abundance domains cov- 
ered by each cluster family and dn/d[Z/H] = N~ 1 dN/d[Z/H], 
being N the projected areal density of each GCs subpopulations at 
a given galactocentric radius, and dn/d[Z/H] comes from the his- 
togram decomposition. Note that the Sn values do not change with 
galactocentric radius and are solely determined by the Z distribu- 
tion parameters of each cluster and associated stellar subpopulation 
(which we also call "blue" and "red" in what follows). 

Both 7 and S were iteratively changed in order to derive a sur- 
face brightness profile that minimises the rms of the residuals when 
confronted with the observed profiles at galactocentric distances 
larger than 120 arcsec. This last value was adopted since both GCS 
exhibit flat spatial density cores that contrast with the peaked shape 
of the galaxy surface brightness. 

These cores can be understood as the result of gravi- 
tational disruption effects that change the original population 
of GCs in the inner regions of galaxies (see, for example, 
Capuzzo-D olcetta & Tessera 1 19991 and references therein) and, 
presumably, become less important for cluster orbits with larger 
perigalactic values. 



5 COLOUR-METALLICITY CALIBRATION 

We present a new colour metallicity relation based on 198 clusters 
that combines revised data for MW GCs and also metallicity data 
obtained for G Cs in three other ga l axies: NGC 3379, NGC 3923, 
and NGC 4649 jPierce et al.l2006ahlPierce et al.l2006bl and Norris 
et al. 2007, in prep.). [Fe/H] values for GCs in these galaxies were 



derived from Lick indices given in the [Thomas et al.l {2004) stellar 
population models. Theses works were selected as they were ho- 
mogeneous both in data handling and in the derivation of the Lick 
indices. 

In the case of MW clusters, we first looked for a trans- 
formation of colours in the Johnson system to (C-Ti). A large 
photometri c samp le, that includes (B-I) colours, is available in 
Reed et alj (fl988l) . In turn, (C-Ti) colours were obtained from 



Harris & Canternal dl977h . Intrinsic colours for these globulars 



were then obtained by using colour ex cesses determin ed by 
iRecio Blanco] d2005h . when available, or fhe lReed et al] Jl988h val- 
ues. As a result we obtain: 



(C - Ti) D = 1.03(±0.02)(B - I) - 0.43(±0.03) 



(11) 



In turn, the extragalactic GCs were observed in the (g-i) Sloan 
colour and transformed to (C-Ti) through two different ways. On 
one side using the (g-i) to (B-I) rela ti on de rived from model in- 
tegrated colours given by iMarastoa d2004h and then to (C-Ti) 
through our own transformation, leading to: 



(C-Tr) 



1.43(±0.03)(g 0.01(±0.02) 



(12) 



Alternatively, Rodge nTet al.l d2006l) have calibrated the Sloan 
indices in terms of Johnson's colour indices that can be transformed 
to (C-Ti) (see, for example. iForbes & Fortel200lh yielding: 



(C - Ti) = 1.39(±0.03)(g -i)+ 0.01(±0.03) 



(13) 



As these transformations are very comparable, within errors, 
we adopted an average of both in order to obtain the extragalactic 
GCs colours on the (C-Ti) scale. 

The intrinsic colours for the extragalactic clusters, were de- 
rived by subtracting the in terstellar reddening excess indicated 
by the ISchlegel et al] J 1998b maps and assuming E(C — Ti) = 
1.97E(B - V). 

The adopted colour-metallicity relation is displayed in Figure 
|7j where a quadratic fit gives a good representation of the data: 



(C - Tx) = 0.94 + 0M8([Fe/H] zn] + 3.5) 2 



(14) 



The no n-linear nature of this relation has been noted by other 
authors (e.g. lHarris & Harrisll2002l : iLee. Lee & Gibson! 1 20021) and 
a linear fit to the data displayed in this last figure leaves significant 
colour residuals both at the high and low metallicity regimes. 

An analysis of colour resid uals as a function of age (for 
globular with ages available in Ide Angeli etafl 120051) or hor- 
izontal branch morphology, t hrough the HB-index given by 
iMackev & van den Be rgh 2005, reveals no trends with these quan- 
tities as displayed in Figure[8] Th ese results are in agr eement with a 
similar analysis presented by Smith & Stradei] d2007l) . who discuss 
those effects for a number of different colour indices. 

We note that the shape of the empirical calibration is simi- 
lar, to within ± 0.015 mags, in (C-Ti), with the colours of the 12 
Gy model, with a Salpeter lumino sity function and blue horizontal 
branch, given in Maraston (2004). This agreement is reached af- 
ter subtracting a zero point difference of 0.065 mags, to their (B-I) 
model colours and then transforming to (C-Ti) through the relation 
given above. Note that these models (shown as triangles in Figure 
are given as a function of total abundance [Z/H]. 

FigureUJalso includes, as a reference, t he NGC 1399 nucleus 
according to the photometry bv lOstrovetSldl993l) and the metal- 
licity {[Fe/H] = 0.4) obtained bv TpicklesI Jl985fK 
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Figure 7. (C-Ti) colour versus metallicity ([Fe/H] on the Zinn and West 
scale) relation derived from Milky Way clusters (filled dots) and "lobulars 
in NGC 3379, NGC 3923 and NGC 4649 (open dot s). The star repre sents 
the nucleus of NGC 1399, triangles are models from Maraston |20Q4|) (see 
text). The continuous line is a quadratic fit adopted as the mean calibration. 
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Figure 8. Upper panel: (C — Ti) colour residuals from the mean colour- 
metallicity calibration for MW clusters included in Figure|7]as a function of 
normalised ages available in lde Angeli et alj (2005;!). Lower panel: (C-Ti)o 
colour residuals for MW clusters incl uded in Figure[7]and horizontal br anch 
morphologies (HB index) available in Ma ckev & van den Berghl42005l) . No 
significant trends are detectable 



Table 3. r 1 / 4 law fits to logarithmic areal densities for globular clusters 
(Ti =21.0 to 23.0) 



colour range 


a (slope) 


b (zero point) 


rms 


NGC 1399 








1.55-2.30 


-0.77 ± 0.03 


3.64 ±0.11 


0.02 


0.90- 1.55 


-0.43 ±0.10 


2.28 ± 0.41 


0.08 


0.90- 1.26 


-0.25 ±0.11 


1.20 ±0.43 


0.08 


NGC 4486 








1.52-2.30 


-0.91 ±0.08 


4.26 ± 0.34 


0.06 


0.90- 1.52 


-0.46 ± 0.04 


2.87 ± 0.20 


0.04 


0.90- 1.21 


-0.22 ± 0.03 


1.64 ±0.12 


0.02 



6 THE GLOBULAR CLUSTERS PROJECTED AREAL 
DENSITY DISTRIBUTIONS 

As shown in FFG05 (figure 5) the slope of the areal density of 
the bluest GCs (i.e., bluer than the blue peak at (C-Ti)=1.26) in 
NGC 1399 is significantly shallower than that corresponding to the 
"whole" blue population (i.e., all clusters bluer than the colour val- 
ley at (C-Ti)=1.55) in the inner region of the galaxy. At larger 
galactocentric radii, these slopes become identical within the un- 
certainties. The significance of that result is analysed in this section 
on the basis of the new data set presented in this work. 

First, we focus on the areal density distributions of clusters in 
the inner regions of both galaxies. The size of this region was de- 
fined aiming at: a) including a large number of cluster candidates; 
b) keeping the overall completeness level of the sample at ~ 95% 
; c) minimising the fraction of contaminating non-resolved field 
interlopers (19% and 11% for NGC 1399 and NGC 4486, respec- 
tively). 

These requirements are met within a circular annulus with in- 
ner and outer radii of 120 and 420 arcsec. Within 120 arcsec the 
searching routines are affected by the galaxy halo brightness while, 
further out in galactocentric radius, the background level increases 
and the effective areal coverage of our images decreases. 

We also set a magnitude range (Ti=21.0 to 23.2) for two 
reasons. On one side, in order to avoid the eventual presence of 
very bright objects whose nature might be connecte d with Omega 
Cen-l ike objects or compact dwarf galaxies (e.g. see lPhillips et al.l 
l200lh . On the other, the GCs colour distribution becomes "uni- 
modal" in NGC 1399 (Ostrov et al. 1998) making difficult a sepa- 
ration between blue and red GCs. 

Given the relatively small angular scale of this analysis we 
adopt r 1//4 laws in order to obtain least squares fits to the logarith- 
mic surface densities within concentric circular annuli (one arcmin 
wide): 

log(den) = ar 1 ^ 4 + b (15) 

The resulting slopes and their associated uncertainties are 
listed in Tableland depicted in Figures [9] The upper two fits, in 
each panel, belong to the red and blue globulars defined in terms 
of the colour valleys at (C-Ti)=1.55 and 1.52 for NGC 1399 and 
NGC 4486, respectively. These fits, that correspond to the regions 
with the highest GCs areal densities, are later overlapped (Figures 
11, 12 , 17 and 18) with profiles that extend to larger galactocentric 
radii. 

In turn, the lower fits belong to GCs bluer than the respective 
blue peaks (at C-Ti=1.26 and 1.21). Both galaxies show very sim- 
ilar behaviours in in the sense that the bluest globulars exhibit sig- 
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nificantly shallower density slopes as found in FFG05 for the case 
of NG C 1399 but using Washington photometry from lDirsch et al.l 
d2003h . 

A further analysis, adopting different colour windows, shows 
no meaningful differences in the slopes of the clusters redder than 
the colour valley, and we consider that they "genuinely" belong to 
a single population. 

The intermediate slope value observed for the whole blue GCs 
sample (compared to the bluest GCs) tentatively suggests that an 
overlap between the blue and red globular subpopulations may oc- 
cur in the colour range defined between the blue peak and the colour 
valley. This overlap would increase the density slope of the so far 
called blue clusters as result of the presence of the blue tail of the 
red subpopulation within their formal domain (i.e., objects bluer 
than the colour valley). 

That effect should decrease with increasing galactocentric ra- 
dius, as the presence of the red clusters becomes less prominent due 
to the steeper density profile of these clusters. This tentative picture 
is discussed in the following sections. 



T 
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Figure 9. Upper panel: Projected areal density for clusters in NGC 1399. 
The upper line belongs to clusters redder than (C-Ti)=1.55 ("red" globu- 
lars). The intermediate line corresponds to clusters with colours between 
0.9 and 1.55 ("blue" globulars according to the most usually adopted defi- 
nition). The lowest line belongs to clusters bluer than the blue peak at (C- 
Tr)=1.26 (or "genuine" blue clusters according to the text). Lower panel: 
Projected areal density for clusters in NGC 4486. The upper line belongs to 
clusters redder than (C-Ti)=1.52 ("red" globulars). The intermediate line 
corresponds to clusters with colours between 0.9 and 1.52 ("blue" globu- 
lars according to the most usually adopted definition). The lowest line be- 
longs to clusters bluer than the blue peak at (C-Ti )=1 .21 (or "genuine" blue 
clusters according to the text). 



7 THE CASE OF NGC 1399 

1) Colour histogram decomposition. 

The background corrected GCs colour histogram is compared 
in Figure \10\ with a synthetic one derived through the modelling 
described in Section|4] This histogram includes ~ 1000 GCs can- 
didates with Ti=21.0 to 23.2 and (C-Ti)=0.90 to 2.30. 

The decomposition process yields 620 clusters to the red sub- 
population with an abundance scale factor Zs(red) = 1.45 ± 
Q.IZq. The remaining 380 globulars are identified as belonging to 
the blue population with an abundance scale Zs(blue) — 0.045 ± 
O.OIZq. Figure [TO] (lower panel), with comparison purposes, also 
displays the Gaussian components that give the best fit to the ob- 
served histograms (blue clusters: (C — 21) = 1.26, at = 0.12; 
red clusters: (C - Ti) = 1.77, oy = 0.20). These fits decrease the 
number of red clusters and increase the number of blue ones sug- 
gesting a smaller degree of colour overlapping between both GCs 
subpopulations in comparison with the results from the model. 

As discussed below, the eventual inclusion of a blue tilt com- 
parable with that adopted for the NGC 4486 blue GCs does not 
have a detectable effect on the shape of the colour histogram. 

Figure QTJJ suggests that a single abundance scale parameter 
for the red GCs population falls somewhat short in the extreme red 
end of the colour histogram. About 5% of that population appears 
definitely redder than the model prediction. A tentative explanation 
might suggest some degree of field contamination in that colour 
range or a possi ble effect con nected with a variation of the [a / Fe] 
ratio with age dKravtsovll2007l) . 

Figure [TITJ also shows that the model colour distribution of the 
red GCs exhibit a "blue" tail (i.e., clusters bluer than the colour 
valley at (C-Ti)=1.55), representing about 31% of the total number 
of red clusters. 

These objects appear as "contaminating" the formal domain 
of the genuine blue globulars and will have an impact on the den- 
sity slopes derived for this last population if only the colour valley 
is adopted as a discriminating criteria between both GCs subpopu- 
lations. 

Alternatively, the model blue GCs barely reach the colour 
valley, suggesting that this population will not affect the estimate 
of the areal density slope of the red clusters if only clusters redder 
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Figure 11. Large scale areal density distribution for globular cluster bluer 
than the blue peak at (C-Ti )=1 .26. This colour domain is practically uncon- 
taminated by the blue tail of the red globulars population and considered as 
the real distributio n of the "genu i ne" bl ue globulars. The filled black dots 
come from data in lBassino et al.l 120061) . The continuous line is projected 
areal density described in text. The short straight line is the fit for the bluest 
globulars in our photometry as depicted in Figure [9] 




Figure 10. Globular populations fit to the background corrected colour his- 
togram for clusters in NGC 1399 is depicted in the upper panel. The vertical 
line at (C-Ti)=1.55 is the colour "valley" usually adopted as the boundary 
between both subpopulations. The lower panel shows each of the globular 
populations derived from exponential distributions in the abundance domain 
Z. Vertical lines indicate the colours peaks of each clusters populations. 
Dotted lines in the lower panel show the components of the best Gaussian 
fit. 



than the colour valley are included in the sampling. 

2) The areal density distribution of the blue and red globulars. 

Due to the features discussed in the previous item, we only 
take clusters bluer than (C-Ti)=1.25 (the blue peak) as tracers of 
the surface density of the genuine blue GCs. Figure [10] shows that 
there would still be a small degree of contamination by the bluest 



clusters of the red population (about 5% of the total sample within 
that colour range). 

The density run with galactocentric radius of these clusters is 
depicted in Figure [TT] In this case we only includ e objects with 
Ti=21 .0 to 23.2 taken from the photometric work bv lBassino et al.l 
(2006), that reaches a galactocentric radius of 40 arcmin. The short 
straight line represents the density fit discussed in Section [6] while 
the continuous line comes from projecting on the sky a volumetric 
density profile: 



p{a) = C(l.0 + (a/rs)Y 



(16) 



where a is measured along the galaxy major axis, a scale length 
rs — 375 arcsec and a spatial cut-off at a galactocentric radius of 
450 kpc. 

The large scale density distribution for the red GCs was then 
derived using only clusters redder than (C-Ti)=1.55 as tracers of 
that population and is shown in Figure[T2] The short straight line is 
the fit discussed in Section|6]while the continuous line is a Hubble 
profile with a core radius r c =60 arcsec. 

In order to estimate the density distribution of the total num- 
ber of clusters for each population, the density of the tracer GCs 
should be increased by factors that take into a account the total 
colour range covered by these populations (adding the blue clus- 
ters redder than the blue peak and the red clusters bluer than the 
colour valley, respectively), as indicated by the colour histogram 
modelling, and also the sampled fraction of globulars within their 
resp ective integrated lumin osity functions. 

iGrillmair et alj |l999) have derived the luminosity functions 
of both blue and red GCs populations in NGC 1399 on the basis 
of HST WFPC2 observations. Assuming fully Gaussian luminosity 
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Figure 12. The same as figure 8 but for clusters redder than (C-Ti)=1.55. 
The continuous line is a Hubble profile with a core radius r c =60 arcsec. 
The short straight line is the fit for the red globulars in our photometry as 
depicted in Figure[9] 

functions, and transforming (B-I) colours to (C-Ti), their results 
lead to turn overs at Ti =23.40 and Ti=23.45 with dispersions of 
1.24 and 1.16 mags for the blue and red populations respectively. 

The combined colour and luminosity completeness factors 
are then, 4.28 for the blue globulars and 3.29 for the red GCs. 

3) The surface brightness profile. 

Surface brightness photometry for NGC 1399 in the B band 
up to a galactocentric radius of 775 arsecs was presented in FFG05 
and compared with other profiles available in the literature (e.g. 
ISchombed[l986l : ICaon et alll999h . 

The predicted blue profile was obtained through the proce- 
dure described in Section [4] and adopting a distance modulus (V- 
Mv) =31.4, corresponding to 19 Mpc (see FFG05 and references 
therein), and an interstellar colo ur excess E(B-R)=0.011 (trans- 
formed from lSchlegel et al.ll 19981) . 

Azimuthal counts do not show a detectable flattening of the 
NGC 1399 GCS and therefore we adopted the average flattening of 
the galaxy, q=0.86, as representative for both cluster populations. 

The model surface brightness profile delivered by ELLIPSE 
is compared with the FFG05 observations in Figure [L31 and corre- 
sponds to 7 = 0.82(±0.05) x 10" 8 and 5=1.1 ± 0.1. The overall 
rms of the fit is ± 0.035 mags. 



8 THE CASE OF NGC 4486 

1) Colour histogram decomposition. 

The GCs background corrected colour histogram is com- 
pared with the model fit in Figure [14] In this case, 800 clus- 
ters were assigned to the red population with an abundance scale 
Zs(red) = 0.90 ± 0.1Z® and 1000 clusters to the blue popu- 
lation with Zs(blue) = 0.012 ± 0.005 practically unde Zq (see 
below). Figure 14 also shows the Gaussian components (blue clus- 



Figure 13. Observed B surface brightness profile for NGC 1399 (FFG05; 
filled dots) confronted with the model fit (open circles). This last profile was 
obtained from a bi-dimensional model image processed with the ELLIPSE 
task within IRAF. Squares and triangles represent the luminosity associ- 
ated with the "blue" and "red" stellar populations. Note that the model fails 
inside 100 arcsec in galactocentric radius where the globular distributions 
shows a flat core. 



ters: (C - Ti) = 1.21, a = 1.12; red clusters: (C - Ti) = 1.72, 
a — 0.20). Here we also adopt an initial abundance Zi — 0.003^© 
for the red clusters. However, as shown in Section|3] the blue GCs 
display an evident tilt that we associate with a change in abundance 
that correlates with the cluster bright ness and, hence, mass (see 
also fig ure 3 in lBrodie & Stra der 2006). In this case, we find that a 
change in initial abundance as a function of brightness: 

AZ = 0.01(23.2 -Ti) (Z e units) (17) 

reproduces the appearance of the blue GCs colour-magnitude di- 
agram. The mean Z for blue GCs with Ti from 21.0 to 21.25 
mags, is 0.0371Zq while for clusters with Ti from 22.95 to 23.2 
is 0.017Zq. These values are consistent with a mass/metallicity 
scaling relation (where M is the clusters mass): 

Z « Af ' 44 (18) 

somewhat smaller than Z ~ jyf° 55 but compar able t o Z w M° 48 
suggested respectively by lHarris et"ai] d2006t) and IStrader et al.l 
(2006). Figure 3 also suggests that the blue GCs tilt, as noted by the 
last authors, is in fact detectable over the whole magnitude range 
brighter than Ti=23.2. 

As mentioned in Section[3] a tilt is not detected in the case of 
the NGC 1399 GCs. In this galaxy, blue GCs exhibit a considerably 
larger Zs(blue) than in NGC 4486 and we suggest that, this larger 
abundance spread makes more difficult the detection of an eventual 
tilt. As an example, Figure [15] displays the colour magnitude 
diagram for the adopted model in the case of NGC 4486, showing 
the blue tilt (left panel). An increase of Zs(blue) from 0.012 to 
0.05 (comparable to that of the blue GCs in NGC 1399) changes 
substantially the appearance of that diagram and makes the blue 



12 Forte et al. 




Figure 15. Model colour diagrams with Zs(blue) = 0.012Zq (left panel) and Zs(blue) = 0.05Zq. Both models include a blue tilt similar to that discussed 
in the text. However, the larger Zs(blue) adopted in the right panel, makes more difficult the detection of the tilt. 



tilt (included in the model) less evident (right panel). 

2) The areal density distribution for blue and red clusters. 

The fact that the NGC 4486 GCS exhibits a noti ceable flatten- 
ing has been pointed out bv lMcLaughlin et al] dl994h . This feature 
is clearly seen in Figure [T6] which shows azimuthal counts within 
a galactocentric circular annulus with inner and outer radii of 120 
and 360 arcsec, performed using our photometry (Ti=21 to 23.2 
and (C-Ti)=0.90 to 2.30). This figure also displays the results from 
a model that includes red GCs with a flattened spatial distribution 
with q=0.80, and blue clusters with q=0.50 (where q=b/a, is the ra- 
tio of the minor to major semi axis). The details of this model, that 
provides a consistent fit to the observations, are discussed below. 

The adopted flattenings come from the assumption that, if 
GCs trace a given stellar population, they should share the same 
flattening. The inner regions of NGC 4486 (where the red stellar 
component dominates the surface brightness) exhibits a q=0.85 (at 
a=120 arcsec) a nd reaches about q= 0.50 at the outermost detectable 
boundaries (see Mih os et al. 2005), where the blue stellar popula- 
tion should become more evident. 

The density run on a large angula r scale was determine d by us- 
ing Suprime camera observations by iTamura et al.l d2006bl) . These 
authors determine areal density in circular annuli on a rectangular 
strip that extends to the east of the galaxy centre. 

We stress that, as those authors use the colour valley at (V- 
I)=1.10 in their photometry as a discriminant between both cluster 
populations, their so called blue clusters will eventually include the 
blue tail of the red population inferred from the colour histogram 
decomposition. We note that Tamura et al. also find that a projected 
NFW profile gives good representation of the areal density distri- 
bution of the so defined blue GCs as FFG05 did in NGC 1399 using 
the same definition for the blue clusters. 

In order to test the compatibility of their observations with 
our approach, we generated a model that assumes that both cluster 
populations follow elliptical distributions with the flattenings men- 
tioned before and a major axis coincident with that of the galaxy 
halo (P.A. ~ 155 degrees). 

The density distribution of the genuine blue clusters in the in- 
ner regions of the galaxy was fit using a surface density profile sim- 



ilar to that adopted for NGC 1399 but with a scale length rs=350 
arcsec. This fit gives an adequate representation to the density de- 
picted in Figure[9](lower panel). 

In turn, model red clusters were generated adopting a lowered 
Hubble d ensity profile (or analytical King profile) with rc=60 arc- 
sec (from|Kundu & Whitm ordl 19981) . In this case, the tidal radius 
was changed iteratively until the best fit to the Tamura et al. densi- 
ties was obtained, yielding rt=3600 arcsec. 

Model GCs colours were generated as described above while 
the Ti magnitudes were derived by adopting fully Gaussian inte- 
grated luminosity functions with turn-overs at Ti=22.9 and 23.2 
and dispersions of 1.38 and 1.55 mags, for the blue a nd red globu- 
lars, res pectively. These parameters were taken from lTamura et al] 
(2006 a), who give values in the V band, and transformed to the Ti 
band (through (V - R) = (V - Ti) = 0.21(C - Ti) + 0.19). 

The completeness factors, that allow an estimate of the total 
number of GCs in each subpopulation from the fractional sampling 
in colour and magnitude, were 3.20 for the blue GCs and 3.31 for 
the red GCs. 

The combined cluster population was then sampled in circular 
annuli, in order to compare with the Tamura et al. density profile 
and taking those GCs bluer than the colour valley at (C-Ti)=1.52 
(or (V-I) « 1.10, following their definition of the blue population). 

The result for this blue population (shifted in log (dens.)) is 
shown in Figure [T71 that also includes a straight line that represents 
the NFW pro file (rs=226 arcsec; 2 kpc at their adopted galaxy 
distance) fit bv lTamura et al. (2006a) to this cluster population. The 
lower line in this diagram corresponds to the genuine blue GCs and 
comes from sampling, also in circular annuli, the projection on the 
sky of an oblate ellipsoid (with q=0.5). This ellipsoid follows the 
blue GCs spatial density dependence mentioned above with a cut 
off at 450 kpc from the galaxy centre. 

Figure[T7]in fact shows that the model discussed in this section 
is able to match the Tamura et al. density profile, which does not 
show a flat core. 

The difference between this last profile and the adopted one 
for the genuine blue GCs, can thus be explained as the result of 
including the blue tail of the red cluster population, characterised 
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Figure 14. A globular population fit to the background corrected colour 
histogram for clusters in NGC 4486 (see Figure |4} is depicted in the up- 
per panel. The vertical line at (C-Ti)=1.52 is the colour "valley" usually 
adopted as the boundary between both populations. The lower panel shows 
each of the globular populations derived from exponential distributions in 
the abundance domain Z. Vertical lines indicate the colour peaks of each 
cluster populations. Dotted lines in the lower panel show the components 
of the best gaussian fit. 



by a steeper spatial distribution, within the sample of GCs bluer 
than (V-I)= 1.10. 

Figure [T8] in turn, shows the comparison of the model with 
the red GCs as defined by Tamura et al. which, not suffering the 
colour overlapping effect, is directly comparable with our model 
red clusters. 



Figure 16. Azimuthal counts (within 22.5 degrees bins) for the globular 
clusters brighter than Ti=23.2 and (C-Ti) from 0.9 to 2.3 (large filled dots). 
The statistical uncertainty of the counts is shown with bars. Open dots come 
from a composite model that includes red clusters with a flattening q=0.8 
and blue clusters with q=0.5 according to text. 
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Figure 17. Projected areal distribution for NGC 4486 GCs bluer than the 
colour valley. Filled dots come from counts within circular annuli given by 
iTamura et al] J2006bl) . Open dots are from the model described in the text. 
The straight line is the NFW fit given by those authors. Note that, following 
this colour definition, the so called "blue" globulars do not exhibit the flat 
inner core. The adopted distribution of the "genuine" blue globulars is also 
shown (lower curve). The dashed line has the slope shown in figure [9] for 
GCs bluer than (C-Ti)=1.21. 



3) The surface brightness profile. 
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Figure 18. Projected areal distribution for NGC 4486 GCs redder than the 
colour valley ("red " ^lobulars). Filled do ts come from counts within circu- 
lar annuli given bv lTamura et al.l i2006bl) . Open dots come from the model 
described in the text. The model counts have been shifted vertically in order 
to take into account the limiting magnitude in that work. The dashed line 
has the slope shown in Figure 9 for GCs redder than C-Ti=1.72. 



Figure 19. Observe d B surface brightness profile for NGC 4486 ( 
ICarter & Dixonll 19781 ; filled dots) confronted with the model fit (open cir- 
cles). This last profile was obtained from a bi-dimensional model image 
processed with the ELLIPSE task within IRAF. Squares and triangles repre- 
sent the luminosity associated with the "blue" and "red" stellar populations. 
Note that the model fails inside 100 arcsec in galactocentric radius where 
the globular distributions shows a flat core. 



Two blue surface brightness profiles with a relatively large 
angular coverag e are a vaila ble for NGC 4486 in the literature: 
ICarter & Dixonl dl978l) and ICaon et alj dl999l) . These profiles, 
along the major axis of the galaxy, show good agreement up to a 
~ 600 arcsec where the Caon et al. profile becomes systematically 
fainter . A comparis on with the diffuse light map in the Virgo clus- 
ter bv lMihos et alj d2005h . in turn, indicates a V surface brightness 
V ~ 26.5 mags, per arcsec'"' at a m 1800 arcsec that imply B= 27.1 
to 27.5 mags, per arcsec'"' which is consistent with the Carter & 
Dixon profile which we adopt in what follows. 

Figure[l9]shows the best fit profile obtained through ELLIPSE 
from the blue synthetic image. The profile corresponds to a dis- 
tance modulus (V-Mv) o =31.0 and a n interstellar reddening E(B- 
V)=0.022 from lSchlegel et alj dl998h - 

The profile fit requires 7 = 1.18(±0.05) x 10~ 8 and 5 = 
1.2 ± 0.1 with and yields an rms of ± 0.07 mags. Again, and as 
already noted for NGC 1399, the flat core of the GCS does not 
allow a proper representation of the inner region of the galaxy. 

The ouput from ELLIPSE shows that, as a result of composing 
two diffuse populations with different Battenings, the galaxy model 
flattening varies with galactocentric radius. This tr end is compared 
with t he ellipticity values (e=l-q) obtained by ICarter & Dixonl 
d 19781) in Figure|20] The overall agreement is acceptable although it 
could be improved if the possibility of a variable q (for one or both 
cluster subpopulations) is allowed. However, the statistical uncer- 
tainties connected with the azimuthal counts prevents a meaningful 
estimate of this eventual dependence. 
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Figure 20. Ellipticity (e = 1 — g) va riation of the NGC 448 6 stellar halo as 
a function of semi-major axis a from lCarter & Dixonl Q978) (open circles), 
compared with the expected variation from the model fit (filled circles). 



9 DEPENDENCE OF RESULTS ON UNCERTAINTIES OF 
THE FITTING PARAMETERS 

The overall results from this modelling in terms of specific frequen- 
cies, characteristic t* parameter (integrated over metallicity for 
each of the cluster populations), diffuse stellar mass and (M /L)b 
ratios are listed in Table @j. 

The total stellar masses given in this table include a correction 
that takes into account the region within 120 arcsec in galactocen- 
tric radius, where the model does not provide an adequate fit. 

In what follows we describe the uncertainties of these results 
in terms of the fitting parameters, 7 and S, as well as those 
connected with the colour-metallicity relation, (M/L) ratios, age, 
and adopted abundance scale. 
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Table 4. Model results for NGC 1399 and NGC 4486. 

NGC 1399 NGC 4486 



Adopted (V-Mv)o 


31.4 


31.0 




Zs(blue) 


0.045 


0.012 


a 


Zs(red) 


1.45 


0.90 




7 


82xKT 8 


1.18xl0~ 8 




<5 


1.10 


1.20 




number of blue globs. 


3900 


7000 


b 


number of red globs. 


4500 


4800 


b 


Sn* (blue globs) 


12.1 


27.9 


c 


Sn* (red globs) 


5.3 


8.4 


c 


t* (blue globs) 3 


44xl0~ 8 


8.41 xlO -8 


d 


t* (red globs) 


85xl0" 8 


1.43X10" 8 


d 


t* (total) 


1.3xl0~ 8 


2.82xl0~ 8 


d 


(M/L) B (blue pop) 


4.2 


3.9 




(M/L) B (red pop) 


9.6 


8.7 




Total stellar mass (Mq ) 


7.2X10 11 


4.8X10 11 


e 


Frac. mass (blue pop.) 


0.18 


0.20 




Frac. mass (red pop.) 


0.82 


0.80 





a) Plus a blue "tilt": AZ = 0.01(23.2 - Ti) 

b) Inside a=1500 arcsec, assuming Gaussian LFs 

c) Intrinsic values defined in terms of their associated stellar 
luminosities in the V band. 

d) Integrated values defined in terms of their associated stellar masses. 

e) Inside a projected galactocentric radius of 100 kpc 



-7 parameter: 

Figure [21] depicts the dependence of both 7 and total pro- 
jected stellar mass (within a constant galactocentric radius of 100 
kpc) with distance modulus. The adopted distance moduli for 
both galaxies are also shown along with a (formal) associated 
uncertainty of ±0.25 mags. Even such small uncertainties, do 
not rule out that both galaxies might be at similar distances 
from the Sun, and in this case, the 7 parameter and total stel- 
lar masses of NGC 1 399 and NGC 4486 would also be very similar. 

-8 parameter: 

8 is independent of the adopted distance and a variation of the 
order of the fit uncertainty (±0.1), mainly impacts on the total mass 
of the diffuse stellar population associated with the blue GCs, that 
changes by ±15%. 

Larger 8 values imply a decrease in the mass of these stars and 
also redder integrated colours of the composite stellar population. 

-M/L ratio and age: 

The (M/L) ratio depends on the age and metallicity of the seed 
GCs. We tentatively adopt an age of 12 Gy comparable to tha t of 
the Milky Way System (see, for examp le.lde Angeli et alj|2005l) . 

Synthetic population models bv lWoirthevI 1 19941) show that a 
variation of ±2 Gy around the adopted model age increases or 
decreases those ratios by 15% without changing the shape of the 
functional dependence with metallicity. Accordingly, masses also 
change in the same proportion. Age variations of that order how- 
ever, do not have a noticeable impact on either the shape of the 
brightness profile or its integrated colour. 

The characteristic (M/L) b ratios for both diffuse populations 
(as well as for the composite stellar population) were obtained by 
integrating over the whole range of abundances determined by their 
respective abundance scale lengths and adopting an upper cut off 
of 4Zq . These results are not critically dependent upon this formal 
upper limit. 

In particular, we note that the (M/L)b values of the red 




Figure 21. Variation of the 7 parameter (lower curves) and total projected 
stellar mass within 100 kpc (upper curves) obtained from model fits as 
a function of the adopted distance modulus (solid line: N1399; dashed: 
NGC 4486). Horizontal lines indicate a change of ±0.25 mags, around the 
adopted distance moduli. 



stellar populations are comparable to that obtained by Sagl ia et all 
d2000h ((M/L) B = 10 ) in the case of the central regions of NGC 
1399 where the red population dominates the integrated luminosity. 

-Abundance scale: 

The relation between [Fe/H]zw and [Z/H], adopted as a con- 
stant over the w hole abundance ran ge, might not be totally appro- 
priate since the Men del et alj d2007l) work does not include GCs at 
a very low abundance regime. We also attempted models adopting 
the Mendel et al. [Z/H]-[Fe/H]zw off set value (0.131) for the red 
GCs and a tentatively larger value (0.3) for the blue clusters. This 
modification leads to larger (M/L)b ratios for the blue popula- 
tion and to an increase of the total mass of about 20% (although 
no significant improvement of the fits of the colour histograms is 
obtained). 



10 IMPLICATIONS OF THE PROFILE FIT 

The 7 and 8 parameters that provide the best fit to the shape of the 
B band surface brightness profile of each galaxy will also lead to a 
given: 

a) Galactocentric colour gradient of the galaxy halo. 

b) Colour off-set between GCs and galaxy halo. 

c) Behaviour of the cumulative GCs specific frequency with 
galactocentric radius. 

d) Metallicity distribution of the diffuse stellar population. 

A comparison of these predicted features with the observed ones, 
as follows, then may provide independent clues about the success 
of the model. 

a) Galactocentric colour gradient of the galaxy halo. 
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Figure 22. (B-R) colour gradient as a function of galactocentric distance for 
the halo (open squares) compared with the mean globular cluster colours 
from the model fit (open circles). The globular mean observed colours, de- 
rived from the photometry presented in this paper, are also shown (f illed cir- 
cles). The short straight line indicates (B-R) colours derived from lMichardl 
l2000l) . Upper panel: NGC 1399; Lower panel: NGC 4486. The dashed lines 
indicate, for each galaxy, the peak colours of the blue and red GCs. 



The main implication of the model is that, at galactocentric 
distances larger than 120 arcsec, the main driver of the galaxy 
colour gradient is the (luminosity weighted) composition of the as- 
sociated blue and red diffuse stellar populations. 

FFG05 presented the expected colour (B-R) gradient for the 
NGC 1399 halo on the basis that the colours of the diffuse stellar 
populations could be identified with the colours of the peaks of 
the associated GCs. That assumption is no longer necessary in this 
work as the colour of each mass element connected with a given 
globular, as well as its mass to luminosity ratio, are determined by 
the metallicity of the "seed" cluster. 



b) Globulars-galaxy halo colour off set. 

The GCs mean integrated colours (including all clusters) in 
massive elliptical galaxies usually exhibit a galactocentric colour 
radient comparable to that of the g alaxy halo but blue ward shifted 
Strom et alll98ll ; lForte et aljl98lh . 

In the context of the model discussed in this work, the cluster 
gradients arise as a consequence of (number weighting) averaging 
the two GCs subpopulations, characterised by different spatial scale 
lengths. The same reasoning apply to the associated diffuse stellar 
populations but, in this case, weighted through the (M/L) ratios de- 
termined by metallicity, then leading to the observed colour off-set. 

Colour gradients derived from the profile fits are shown in 
Figure [22] In these diagrams, the predicted halo colours are com- 
pared with the mean globular integrated model colours and also 
with those obtained from the photometry presented in section II. 
These diagrams show that, in fact, the galaxy halos exhibit colour 
gradients comparable, but redder, than those of the GCs. 

Figure [22] als o show s that the (B-R) colour gradients deter- 
mined by MichanJ feOOd) in the inner regions of the galaxies are in 




Figure 23. Cumulative globular cluster specific frequencies for NGC 1399 
(lower curve) and NGC 4486 (upper curve) derived from the model fits. The 
vertical lines indicate galactocentric radius of 40 kpc (dotted line: NGC 
1399; dashed line: NGC 4486) and indicate Sn fa 3.5 and Sn fa 8.5, 
respectively. 



very good agreement with with the predicted colours of the halo. 

c) Cumulative globular cluster specific frequency. 

Figure [23] shows the galactocentric variation of the cumulative 
specific frequency derived from the best fit models. The GCs popu- 
lations inside a gala ctoc entric radius of 120 arcsec w ere taken from 
iForbes et alj i 19981) and lKundu & Whitmord dl998h . Even though 
each cluster subpopulation has its own intrinsic frequency, a varia- 
tion of the composite Sn is expected as the number ratio of blue to 
red GCs changes with galactoce ntric radius. 

The parametric Sn values ( McLaughlin et a i] |l994t) . defined 
at a galactocentric radius of 40 kpc, from this figure are Sn ~ 3.5 
and Sn fa 8.5 for NGC 1399 and NGC 4486, respectively. These 
values are considerably low er than previous es timates given in the 
literature as already noted in iForte etaljj2002l) . 

d) The metallicity distribution of the diffuse stellar population. 

The shape of the [Fe/H] distribution expected for the diffuse 
stellar population was schematically derived in FFG05 on the basis 
of estimating a characteristic (M/L)b ratio and intrinsic specific 
frequency for each cluster population. In contrast, in this work each 
stellar mass element has a given metallicity, and hence (M/L) ratio. 
The statistic distribution of these masses, as a function of [Fe/H] is 
given in Figure 24 for NGC 1399 and NGC 4486. For both galax- 
ies we show the inferred metallicity distribution, for galactocentric 
ranges from 10 to 15 and 15 to 25 kpc convolved with a Gaus- 
sian kernel (dispersion: 0.20 dex) aimed introducing some degree 
of smoothing comparable to observational errors. 

These mass statistics can be transformed into star-number 
ones under the assumption that the stellar luminosity functions do 
not depend stron gly on metallicity. A c omparison with the cases 
of NGC 5 128 dHarris & Harris! [2001 or lReikuba etai]|2005h and 
also M31 ( IDurrell et alj|200ll) shows good qualitative agreement, 
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Figure 24. Stellar mass histogram as a function of metallicity for NGC 
1399 (upper panel) and NGC 4486 (lower panel). These histograms belong 
to elliptical galactocentric radii between 10 and 15 kpc (solid line) and 15 
to 25 kpc (dotted line) at an adopted distance of 19 Mpc (NGC 1399) and 
15.8 Mpc (NGC 4486). The histograms asume an a ratio of 0.3 



i.e, the presence of a broad high metallicity component and an ex- 
tended tail towards low metallicity that become s more evident as 
galactocentric radius increases. More recently, Mou hcinel J2006I) 
presents stellar number statistics for a number of edge on spirals 
which also shows low metallicity skewed distributions, a feature 
that seems independent of the galaxy morphology. 



11 CAVEATS ABOUT THE MODEL 

The model described in this work has several caveats, namely: 

-colour bimodality is attributed to two different GCs subpopu- 
lations. An alternative view, based on the presence of an inflection 
region in the colour-metallicity relation as the main d river of the 
shape of the colour histograms, has been suggested bv lYoon et al.l 
(2006). So far, however, neit her our calibration n or recent spectro- 
scopic results in NGC 4472 JStrader et alj|2007h seem to support 
th at situation. Fur t her re sults in this last direction can also be found 
in lKundu & Zepl J2007t) for the case of NGC 4486. 

- Although an exponential dependence of the number of GCs 
with abundance was adopted, more complex functional depen- 
dences, hidden by the noise of the GCs statistics, could not be ruled 
out. 

- An abundance dependence with brightness, that leads to a 
blue tilt in the case of the blue GCs, is included in the models. 
However, we cannot infer whether the tilt is just a local effect or if 
it is actually shared by the diffuse stellar population associated with 
the blue clusters. Nevertheless removing the tilt from the models, 
has little impact on the output integrated colours that then would 
become slightly bluer (w 0.015 mags, in (B-R)). 

-The two dominant cluster subpopulations are assumed to be 
coeval. Further refinement of this aspect could be incorporated 
onc e meaningful ages become available. On one side some works 
(e.g I Jordan et alj 2002) do not find detectable age differences be- 
twee n the blue and red populations in NGC 4486. On the oth er, 
(e.g. [F orbes et al. 200ll |pierce et al.ll2006al : [pierce et alj|2006bl or 
iHempel et alj|2007h do find a certain fraction of intermediate age 
clusters in NGC 1399 and in other galaxies. 

-The [a/Fe] is adopted as constant and equal for the clus- 
ter subpopulations in both galaxies in order to derive the [Fe/H] 
stellar distributions . Even though t he adopted value i s appropri- 
ate for the MW (see lCarnevl l 19961 orlThomas et alj|2003h and also 
representative for ellipticals jPuzia et alj|2005r) . a possible varia- 
tion with metallicity as noted by these last authors and, earlier, by 
IShetroneTC ote & Sargent] J200lh . cannot be ruled out. 



12 DISCUSSSION AND CONCLUSIONS 

The results presented in previous sections show that the surface 
brightness profiles of both NGC 1399 and NGC 4486 can be traced 
using a common link between GCs and the stellar halo popula- 
tions in these galaxies. This link imply that the number of GCs per 
diffuse stellar mass, defined as t = 7 exp (—5 [ Zl H] ) .increases 
when chemical abundance decreases. We note that lHarris & Harrisl 
d2002h had already found an increase of Sn with decreasing GCs 
metallicity in the case of NGC 5128. 

This suggests that, on a large scale, the dominant globular 
cluster subpopulations formed along major star forming episodes 
and following a similar pattern. However, it is not yet clear whether 
abundance, through the role that it plays in the t parameter, is the 
physical reason that governs the fraction of clustered to diffuse stel- 
lar mass or, eventually, some other "hidden" variable in turn corre- 
lated with abundance. 

The quality of the profile fits is comparable to any other para- 
metric approximation in a range that covers a galactocentric ra- 
dius from 10 to 100 kpc. In the inner regions, GCs fail to map 
the brightness distribution probably as a consequence of cluster de- 
struction processes due to gravitational effects. However, it seems 
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that survivor clusters, with large perigalacticon orbits, are still able 
to trace their associated stellar populations. It is wort h mention- 
ing th at, based on far UV observations of NGC 1399 . iLotz et all 
(2000), find arguments that support the coexistence of two widely 
different stellar populations, in terms of chemical abundance, in the 
nucleus of the galaxy. 

Even though the approach only aims at reproducing the bright- 
ness profiles, other connected features as galactocentric colour 
gradients, GCs-halo colour offset, cumulative cluster specific fre- 
quency and inferred stellar metallicity distributions compare very 
well with observations. 

It seems also remarkable that the common quantitative GCs- 
stellar halo link works in both galaxies although their cluster popu- 
lations exhibit detectable differences in cluster numbers and chem- 
ical abundance. Both GCs subpopulations have larger abundance 
scale lengths (as defined in Section 4) in NGC 1399 than in NGC 
4486 leading to total mean abundance ratios of 1.4 and 2.5 for the 
red and blue GCs, respectively. 

As shown in Figures 10 and 14, the approach presented in this 
paper delivers GCs colour distributions that are not strongly dif- 
ferent from a two Gaussians fit requiring five free parameters, a 
common procedure in the literature (e.g Ostrov, Forte & Geisler 
1998). However, we note that our models indicate a lower ratio of 
the number of blue to red GCs and only require three free parame- 
ters. 

Blue GCs in NGC 4486 have a small abundance scale lenght, 
about four times smaller than that of the blue GCs in NGC 1399, 
and we speculate that this may be connected with a shorter forma- 
tion time scale. In turn, that lower abundance spread may be the 
reason behind the presence of a blue tilt, connected with cluster 
mass, in the colour magnitude diagram. This feature seems absent, 
or probably masked, by the larger abundance scale of the blue GCs 
in NGC 1399, a situation that may also hold in other galaxies (e.g. 
NGC 4472, Strader et al. 2006). 

This last result argues in favour of the idea that blue 
GCs "know" about the gal axy they are associated with 
( Strad er. Brodie & Forbes! [2004). In any case, and in both galax- 
ies, the blue GCs barely reach an abundance close to [Z/H]=-0.5. 
The reason for this upper metallicity cut off may be connected with 
some kind of sincrhonising event as the re-ionis ation of the Uni- 
verse JCertl200lklSantol2003l ; rRhode & Zepj|2005h . However, the 
different abundance scale lengths of the blue GCs also suggest that 
a local phenomenon, distinct for each galaxy, may have also played 
a role in modulating the the star formation rate (e.g. the onset of 
galaxy nuclear activity). 

The overall picture seems co nsistent with some scenarios al- 
ready discussed in the literature jForbesetal]|l997l :l Harri s et al.l 
2006) that invoke two different cluster formation mechanisms and, 
probably, environmental conditions. 

In contrast with the relative abundance homogeneity and large 
spatial (half density) core radii (~ 25 kpc) of the blue GCs, the 
red ones exhibit a large abundance heterogeneity and much smaller 
core radii (~ 5 kpc) also shared by the red diffuse stellar popula- 
tion. This degree of heteregeneity may be connected , for exam ple, 
to mergers of different nature (e.g. lSchweizer & S eitzer 1993]). 

The total globular cluster formation efficiencies, in terms of 
stellar mass, indicated by the models, and adopting an average clus- 
ter mass of 2.5 x 1O 5 M , are 2.4 x 10" 3 for NGC 1399 and 
5.6 x 10~ 3 for NGC 4486. T hey are comparable to the efficiency 
derived bv lMcLaughlinl l l 19991) although, in that case, the definition 
of efficiency included total baryonic mass. 

Blue clusters show a higher formation efficiency (in terms of 



the stellar mass they are associated with), when compared with the 
red ones, probably as a consequence of a lower star formation effi- 
ciency during the early phases of galaxy formation at a low metal- 
licity regime. 

Although there is no strong evidence of an age d ifference be- 
tween the blue and red globular subpopulations (e.g. Ijordan et al.l 
2002) within the uncertainties of the measurements, a possible 
temporal sequence that assumes the formation of the blue popu- 
lation first, cannot be discarded as a result of the relatively small 
time scales involve d at the early phases of galaxy formation (see 
iBeaslev et al1l2003h . In this frame, the chemical enrichment pro- 
vided by a presumably progenitor blue population might be impor- 
tant to boost later stellar formation efficiency through an abundance 
enrichment that may reach [Z/H] ~-0.60 within 100 kpc of the 
galaxy nucleus. 

Some scenarios sug gest that blue GCs formation is asso- 
ciated with dark matt er dBeaslev et al.l |2002| ; iMoore et alj 120061 : 
Prieto & Gnedin 2006), and it is tempting to look for such a con- 
nection. For example, the results listed in Table [4] show that while 
both galaxies have a similar total number of red GCs, NGC 4486 
outnumbers NGC 1399 in a factor of about 1.8 in terms of blue 
GCs. Dark mass estimates within a galactocentric radius of 100 
kpc are 3.4xl0 12 M (T) for NGC 1399 (extrapolating data from 
Richtler et alj|2004l) and 7.4x 10 12 A/ Q for NGC 4486 dCote et al.l 
Il998h . leading to a ratio ~ 2.0 comparable to that in the number of 
blue GCs. 

The similarity of the stellar galaxy masses, and the differ- 
ence i n their total masses, had already been noticed by I Jones et al.l 
dl997h on the basis of their X ray analysis. 

FFG05 found that, adopting their definition of blue clusters, 
the density pr ofiles of the NGC 1 399 GCs could be fit with a 
NFW profile dNavarro et al] 1 19961) with a scale length of 375 
arc sec, coincident wi t h that derived for the inferr ed dark matter 
halo bv lRichtleret alj j2004l) . lTarnura et all d2006al) also perform a 
NFW profile fit to the blue clusters in NGC 4486. However, both 
approaches deserve a revision since, on the basis of the results 
presented in this work, the "genuine" blue GCs exhibit a rather 
extended inner core in their surface density profiles. As shown 
here, these cores have been disguised by the inclusion of the blue 
tail of the red subpopulation within the "blue" GCs sample. This 
overlapping should be even more severe when using colour indices 
less sensitive to metallicity than (C-Ti) and should be taken into 
account when doing, for example, kinematic analysis of the cluster 
subpopulations. 
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